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ABSTRACT 



Context. Clusters of galaxies are believed to be capable to accelerate protons at accretion shocks to energies exceeding lO'^ eV. 
At these energies, the losses caused by interactions of cosmic rays with photons of the Cosmic Microwave Background Radiation 
(CMBR) become effective and determine the maximum energy of protons and the shape of the energy spectrum in the cutoff region. 
Aims. The aim of this work is the study of the formation of the energy spectrum of accelerated protons at accretion shocks of galaxy 
clusters and of the characteristics of their broad band emission. 

Methods. The proton energy distribution is calculated self-consistently via a time-dependent numerical treatment of the shock accel- 
eration process which takes into account the proton energy losses due to interactions with the CMBR. 

Results. We calculate the energy distribution of accelerated protons, as well as the flux of broad-band emission produced by sec- 
ondary electrons and positrons via synchrotron and inverse Compton scattering processes. We find that the downstream and upstream 
regions contribute almost at the same level to the emission. For the typical parameters characterising galaxy clusters, the synchrotron 
and IC peaks in the spectral energy distributions appear at comparable flux levels. 

Conclusions. For an efficient acceleration, the expected emission components in the X-ray and gamma-ray band are close to the 
detection threshold of current generation instruments, and will be possibly detected with the future generation of detectors. 

Key words, shock acceleration - radiation mechanisms: non-thermal 



1. Introduction 

Rich clusters of galaxies are the largest virialised structures in 
the Universe, with ty pical sizes o f a few Mpc and masses up to 
IO'^Mq or more (see lSarazinlll98 8 for a review). In the standard 
picture of cosmic structure formation, the structure's growth 
is driven by gravitational instability. This process is hierarchi- 
cal, with larger systems forming later via the assembly of pre- 
existing smaller structures. Within this scenario, galaxy clusters 
form via mergers, a nd their age can be estimated to be of the 
order of 10 Gyr (e.g. lBorgani & Guzzoll2001i) . In addition, cold 
material from the surrounding environment is continuously in- 
falling, due to gravitational attraction, and an expanding shock 
wave, called the accretion shock, is expected to form at the clus- 
ter b oundary and to car ry outward the information of virialisa- 
tion dBertschingeij 19851) . Numerical simulations have confirmed 
the appearance of the so-called a ccretion shocks during structure 
formation (e.g. lKang et al.lll994l) . 

The detection of a tenuous and diffuse synchrotron radio 
emission from about on e third of rich clusters of galaxies (e.g. 
iGovoni & Ferettill2004l) reveals the presence of a diffuse mag- 
netic field and a population of high energy electrons. Using ex- 
clusively synchrotron data, though, it is possible to get infor- 
mation only on the product of particles density and magnetic 
field energy density (unless adopting assumptions like equipar- 
tition). Detection of non-thermal X-rays have also been claimed 



from a few of such sources and generally interpreted as inverse 
Compton emission from the same populat ion of electrons (e.g. 
iFusco-Femiano et a l 1999: Eckert at al . 2008). Since IC depends 
on the electron density but not on the magnetic field, combining 
the data allows us to break the degeneracy and thus determine 
values of B of the order of 0. l/iG. Significantly higher values, of 
the order of a few microG auss, are obtained from Faraday rota - 
tion measurements dClarke et al. 1 120011: ICarilU & Tavloiil2002l) . 
Thus the uncertainty in the determination of the magnetic field 
is quite large, allowing value s to vary of about one order of mag- 
nitude dNewman et al.ll2002h . 



According to theoretical models, the electrons re- 
sponsible for the radio and non-thermal X-ray emission, 

are produced through diffe r ent accelerati on mechanisms 

iTribbld fT99l [Sarazin 199^ 

20011: iPetrosianI 1200 ll) . Alternatively, the radio 



synchrotron radiation can be emitted by secondary electrons 
produced at in teractions of accelerated protons with the in 
tracluster gas dDennisonl 119 80': Blasi & Colafrancesco '199' 



* work performed as IMPRS fellow at the Max-Planck-Institut fiir 
Kernphysik, Heidelberg - e-mail: giulia. varmoniOcea. fr 



Atovan & V61^ |2000t |Ddag & EnBling 2000; Kushniretal 
20091) . In this scenario, gamma ray emission is expected due 
to the decay of n eutral pions produced in p-p interactions 
dVolk et al.l 1 19961: iBerezinskv etakl Il997l) . The non-thermal 
X-rays emission also can be related to synchrotron radiation 
of secondary electrons, but of much higher energies, produced 
in photon-photon (Timokhin et al. 2004) and proton-photon 
dAharonian II2002I: ilnoue et al.ii2005k) interactions. 



2 



G. Vannoni et al.: Acceleration and radiation of ultra-high energy protons in galaxy clusters 



Several particle acceleration mechanisms have been pro- 
posed to operate in clusters of galaxies (see i Blasi et aTl |2007| 
for a review). In particular, it has been argued that large scale 
shocks can effectively a ccelerate electrons and protons up to 
ultrar elativistic energies ('Norman et aL 1995 : Berezinsky et alj 
19971 IT oeb & Waxman 2000; Mini ati et al.1 120011; iBlasil 120011; 
Gabici&Blasi 2003; Rvu et alj |2003|; iBerrington & Dermerl 
2003; Pfrommer et al. 2006). In a recent paper we performed de- 
tailed calculatio ns of diffusive shock acceleration of electrons in 
galaxy clusters dVannoni et al.ll2069l) . Electrons can be acceler- 
ated up to 100 TeV at cluster accretion shocks, with synchrotron 
X-ray and IC gamma-ray radiation components produced mainly 
in the downstream region. While the maximum energy of elec- 
trons is limited by synchrotron and IC energy losses, the protons 
can be accelerated to much higher energies. 

Indeed, according to the Hillas criterion (iHillasI 11984 . 
galaxy clusters are amongst the few source populations capable, 
as long as this concerns the dimensions of the structure and the 
value of the magnetic field, to accelerate protons up to 10^'^ eV. 
Moreover, clusters are cosmological structures and their lifetime 
is comparable to the age of the Universe, therefore, if acceler- 
ation takes place in such objects, it can continue up to ~ 10'" 
yr. On the other hand, high energy protons up to energies of 
the order of 10'^ eV and possibly higher are well confined in 
the volume of the cluster over this time scale ( !Volk et al.l[T996t 
iBerezin sky et al.l [T997). This results in an effective accumula- 
tion of high energy particles in the cluster. The thermal energy 
budget of rich clusters, estimated from measurements of their 
thermal X-ray emission, is of the order of lO*"^ - 10*"^ erg. The 
non-thermal component seems to be const rained to a few per- 
cents of the thermal energy ~ 10^^ erg (iPerkins et aP l2006t 
lAharonian et"ani2009t lAleksic et al]l2009h . This value is also 
comparable with the magnetic field energy, assuming 1 fiG and 
a spherical cluster of 3 Mpc of radius (e.g. the size of the Coma 
cluster). 

Although the dimensions of the system, the strength of the 
magnetic field, and the age of the accelerator formally allow 
protons to be accelerated up to 10^" eV, the particles lose, in 
fact, their energy via pair and pion production in the interac- 
tions with the photons of the Cosmic Mi crowave Background 
(C MBR) radiat i on fie ld. As discussed in Norman et aljll995h 
and lKang et al.l (Il997l) . for a shock velocity of a few thousands 
km/s and a magnetic field of the order of lyuG, the shock ac- 
celeration rate is compensated by the energy loss rate at ener- 
gies around 10'"^ eV (the exact value depends on the assumed 
dif fusion coefficient and also on the shock geometry as shown 
bv lOstrowski & Siemieniec-Ozie blo 2002). The interactions of 
ultra-high energy protons with the CMBR lead to production of 
electrons in the energy domain (> 10'^ eV) which is not accessi- 
ble through any direct acceleration mechanism. These electrons 
cool via synchrotron radiation and IC scattering on very short 
timescales (compared to both the age of the source and interac- 
tion timescales of protons). For the same reason they are rather 
localised in space, i.e. are "burned" not far from the sites of their 
production. Therefore the corresponding radiation components 
in the X-ray and gamma-ray energy bands are precise tracers of 
primary protons, containing information about the acceleration 
and propagation of t heir "grandparents" . This interesting feature 
has been indicated in lAharonian I (12002!) in the general context of 
acceleration and propagation of ultrahigh energy protons in large 
scale extragalactic structures. More specifically, this issue was 
discussed by Inoue et al. (2005) for objects like the Coma galaxy 
cluster However in these studies an a priori spectrum of pro- 
tons has been assumed in the "standard form" E^^ exp{-E /E*). 



In fact, the energy distribution of protons accelerated in galaxy 
clusters via Diffusive Shock Acceleration (DS A) mechanism can 
be quite different from this simple form, as it is demonstrated be- 
low. 

In this paper we study the process of proton acceleration 
by accretion shocks in galaxy clusters taking into account self- 
consistently the energy loss channels of protons related to their 
interactions with the CMBR photons. We m ake use of the nu- 
merical approach presented in I Vannoni et al.l (^009). In Section 
|2]the calculation is introduced and the accelerated proton spec- 
tra are derived. In particular, we show that electron/positron pair 
production is the dominant energy loss channel for protons. In 
Section[3]we calculate the spectra of secondary pairs. The broad- 
band emission produced by the secondary electrons during their 
interactions with the background magnetic and radiation fields 
is presented in Section |4] In Section |5] we study, in a simpli- 
fied fashion, the impact of shock modification by efficiently ac- 
celerated protons on the acceleration and emission features. We 
briefly discuss and summarise the main results in Section |6l 

2. Proton Acceleration and Energy losses 

The accretion of cold external material onto a hot rich cluster 
of galaxies can lead to the formation of a strong shock at the 
position of the virial radius of the cluster The shock velocity 
can be estimated via the free fall velocity of the infaUing matter 
crossing the shock surface, namely: 
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It is interesting to note that the shock velocities for this case 
are comparable to the shock speeds typical for young supernova 
remnants. 

We approximate the spherical shock locally as plane and 
work in the reference frame where the shock is at rest and the 
plasma moves along the jc-axis, perpendicular to the shock sur- 
face, from -oo far upstream, to +oo far downstream and the 
shock is located at jc = 0. In the following, all the quantities in 
the upstream region will be indicated with the subscript 1, and 
all the quantities downstream with the subscript 2. To obtain the 
accelerated proton spectrum under the effect of energy losses in- 
duced by the interact ion with the CMBR, w e adopt the numerical 
scheme described in I Vannoni et al.l ( l2009l) . The time dependent 
transport equation for the particle distribution function in phase 
space is: 

df{x,p,t) df(x,p,t) d I df{x,p,t)\ pdudf(x,p,t) 

7~ \D{x,p) — 

ot ox ox\ ox ) 3 ox op 
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(p^L(x,p)f(x,p,t))^Q(x,p), 



(2) 



where L(x,p) - -p is the energy loss rate and Q{x,p) is the 
injection term; u represents the bulk velocity of the plasma in 
the shock rest-frame. The momentum and space dependence of 
the energy loss term as well as of the diffusion coefficient can 
be chosen of any form so that the formulation is general. Eq. 
^ is solved numerically, inserting the boundary conditions at 
the shock location and at upstream/downstream infinity. For the 
boundary at the shock we consider again Eq. (|2]l and integrate it 
between x = 0_ immediately upstream and x - 0+ immediately 
downstream, leading to: 
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Fig. 1. Acceleration and energy loss time scales as a function of the 
proton energy. The acceleration time scales are obtained for the values 
of the upstream magnetic field Sj reported in figure and a downstream 
magnetic field Bi= ABi. The thick lines correspond to a shock velocity 
of 2000 km/s, the thin lines to a velocity of 3000 km/s. As an horizontal 
dotted line we report the estimated age of the Universe, for comparison. 



At +00 we set fix, p) - 0. 

We assume that injection happens at the shock surface and 
that it can be described by a delta-function in momentum; 

Q(x,p) = Qq6(x)6(p - pq), 

where po is the injection momentum and Qo a normalisation con- 
stant. 

As pointed out in iNorman et al.l (Il995h and iKang et aTl 
(Il997h . the relevant energy loss channel for the protons accel- 
erated at cluster accretion shocks is the interaction with the 
CMBR. Other radiation fields do not play any significant role in 
this regard. The interaction proceeds via two main processes of 
electron/positron pair production, also known as Bethe-Heitler 
pair production, and photomeson production. In Fig. [T] we show 
the cooling time-scales of protons due to these two processes 
(solid curve). Pair production dominates up to energies of about 
5 X 10'^ eV. Once the particle energy has passed the threshold 
for meson production, this latter becomes rapidly dominant. The 
thick lines in Fig. [1] coiTespond to the acceleration times for a 
shock velocity of 2000 km/s and an upstream magnetic field of 
0.1 juG (dotted), 0.3 fxG (dashed) and 1 fiG (dash-dotted). For the 
same magnetic field, but for the shock speed of 3000 km/s, the 
acceleration times are shown with thin lines. The downstream 
magnetic field is assumed of the form B2 = ^Bi, where ^ is the 
compression factor at the shock. For a strong linear shock, the 
parameter ^ varies between 1 to 4, depending on the orientation 
of the upstream magnetic field. We choose ^ = 4, appropriate 
for Alfvenic turbulence. Fig. [T] shows that if protons can be ac- 
celerated beyond few times 10'^ eV, the cut-off energy as well 
as the shape of the particle spectrum in the cut-off region are de- 
termined by pair production. Note that even for extreme cluster 
parameters, the protons do not achieve energies for which the 
energy losses would be dominated by photomeson production. 

In this work we choose the following values of model pa- 
rameters: a shock velocity of 2000 km/s and a magnetic field up- 
stream Z?i = 0.3 juG, with B2 - 4Z?i (thick solid line in Fig.[T|. 
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Fig. 2. Proton spectra at the shock location for an acceleration time of 
10 Gyr (solid line) and 5 Gyr (dashed) for a shock velocity of 2000 
km/s, a magnetic field upstream Sj = 0.3 /iG and a magnetic field 
downstream 61 = 4B[. 



We do not take into account the evolution of the CMBR with red- 
shift. This simplification is justified by the fact that the acceler- 
ation rate and the energy loss rate become comparable on a time 
scale of about 5 Gyr. Therefore, for local clusters, this time scale 
places the start of the acceleration process at an epoch where 
the redshift is z «: 1 and we can neglect t he evolutio n of th e 
CMB radiation. Moreover, as discussed in iKang et aP (Il996l) . 
the higher temperature and energy density of the CMBR at large 
redshifts inhibits the acceleration of protons to ultrahigh energies 
at early epoch. 

In Fig.|2]we plot the calculated proton spectra at the shock 
location. Our results are more conveniently shown as a function 
of the particle energy, rather than their momentum. We perform 
our calculation for an age of the system of 5 Gyr and 10 Gyr 
The difference between the spectra at the two ages demonstrates 
that the system does not reach steady state. Therefore, a time 
dependent calculation is essential in order to model correctly the 
particle spectrum. The cut-off energy is located around 7 x 10'^ 
eV. A small bump is present in the spectra around the cut-off 
energy, due to the flattening of the energy loss time scale in that 
energy range. The prominence of the feature increases with time. 

Once the particle distribution function is obtained for ev- 
ery point of the phase space, we proceed to calculate the spec- 
trum integrated over space F{E) - J f{x, E)dx, in order to ob- 
tain the production spectra of secondary pairs. The results are 
plotted in Fig. [3] for the up- and down-stream regions of the 
shock separately. The spectral features are evident in the spa- 
tially integrated spectra. At energies below the cut-off, the up- 
stream particle spectrum is hard. This is due to the fact that par- 
ticles can propagate ahead of the shock front over a distance of 
the order of their diffusion length, defined as xo = D{E)/ui. 
For Bohm-type diffusion xo oc E, and the integrated spectrum 
upstream is thus proportional to the particle spectrum at the 
shock multiplied by E. Downstream, on the other hand, at en- 
ergies low enough that energy losses are negligible, the spec- 
trum is f2{x,E) = fo(E) oc E^^. In this component, a bump 
forms around 10'** eV due to the fact that around that energy 
(close to the threshold of the process) the pair production loss 
rate changes behaviour from a very steep dependence on the 
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Fig. 3. Spatially integrated spectra for the proton distributions in Fig. [2] 
for an accelerator age of 10 Gyr (solid lines) and 5 Gyr (dashed lines). 
The lines represent: 1 and 3 the downstream contributions, 2 and 4 the 
upstream ones, and the sum of the two. 



particle energy to an almost flat distribution. Therefore, parti- 
cles at slightly lower energies are only marginally affected by 
losses, as we can see from the fact that, below 10'^ eV, a power 
law spectrum oc E^^ is recovered, as expected for an uncooled 
spectrum; on the other hand, particles at slightly higher energies 
experience a very fast cooling so that there is an effect of ac- 
cumulation arou nd the threshold. The effect is a nalogous to the 
one discussed in lBerezinskv & Grigoreval(ll988l) . explaining the 
appearance of a bump in the ultrahigh energy Cosmic Ray spec- 
trum due to photomeson production during propagation. Quite 
interestingly, as a consequence of the shallow dependence on 
energy of the loss time scale at the cut-off energy and above it, 
the spectrum in the cut-off region is smoother than a simple ex- 
ponential behaviour At energies around 10^" eV, a second steep- 
ening is in fact present (not shown in figure), due to the effect of 
pion production, however at that energy the flux is suppressed 
by several orders of magnitude with respect to the cut-off value 
and thus meson production is negligible for changing the proton 
spectrum. 

It is worth noting that, even though the contribution from the 
downstream region is dominant at low energies, around the cut- 
off the contribution of the two components differ only by a factor 
of few. 
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Fig. 4. Production spectra of pairs. Solid line: downstream spectrum 
for an accelerator age of 5 Gyr; dashed line: upstream spectrum. 
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Fig. 5. Energy loss rates for pairs due to synchrotron and IC processes. 
Solid line: downstream (magnetic field 82= 1.2 //G), dashed line: up- 
stream (magnetic field Sj = 0.3 /iG). 



3. Spectra of Secondary Electrons 

In the following we present the results calculated for a cluster 
with an age of 5 Gyr. For simplicity we assume that the magnetic 
field is homoge neously distributed over the cluster. The formal- 
ism proposed bv' Kelner & AharonianI (12008 ') was used for calcu- 
lations of the production spectra of secondary electron-positron 
pairs. Fig. |4] shows the spectra of secondary electrons up and 
downstream produced by the proton spectra shown in Fig. |3] 

Once produced, the electrons (hereafter the term "electrons" 
implies both electrons and positrons), lose energy at interac- 
tions with the CMBR and the ambient magnetic field. The en- 
ergy loss rates for these two processes are shown in Fig.|5] The 
curves are multiplied by E^^ in order to indicate the transition of 



the Compton losses from the Thomson to Klein-Nishina regime 
which takes place around 10 TeV. The curves correspond to the 
sum of synchrotron and IC energy loss rates in the downstream 
and upstream regions (solid and dashed line, respectively). In 
both cases the target photon field for IC is the CMBR with 
temperature 2.7 K. The magnetic field downstream is 4 times 
larger than upstream. While at low energies the Compton losses 
dominate over the synchrotron losses, at very high energies, 
E > 1000 TeV, because of the Kiein-Nishina effect, the syn- 
chrotron losses become the dominant channel of cooling of elec- 
trons. 
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Fig. 6. Electron spectra after taking into account energy losses. The 
solid and dashed lines refer to the downstream and upstream regions 
respectively and the spectra are normalised to a total energy in acceler- 
ated protons of 10*^ erg in the volume of the cluster, a) Spectra obtained 
considering pure synchrotron cooling in a magnetic field Bi =0.3 fjG 
upstream and St = 1.2 yuG downstream, b) Same as in panel a but with 
the realistic combination of IC and synchrotron losses as in Fig. [5] 



For the given injection (source) rate Qmj{E), the energy 
distribution of electrons F{E) is described, in the continuous 
energy-loss approximation, by the equation: 



dF(E) 
dt 



d IdE \ 



(4) 



Here we assume that the particles do not undergo any fur- 
ther acceleration during their lifetime. In the steady state limit, 
dF(E)/dt - 0, the energy distribution of electrons is given by: 



(5) 



The steady state solution coiTectly describes the spectrum of 
electrons, given that the lifetime of high energy electrons is 
shorter than the age of the source. 

The two panels in Fig. |6] show the cooled spectra of elec- 
trons downstream (solid) and upstream (dashed). The normali- 
sation is obtained assuming a total energy in accelerated protons 




103' ^ 

iiiiJ J llllJ J llllJ J llllJ J I 



102 104 106 108 1010 1012 1014 1016 

£ [eV] 



ra 

~g 10- '3 

o 

lio-'5 



10- 





MMI^ MM^ MMI^ 1 


111^ MMI^ MM^ MMI^ MM^ 


MMI^ MM^ MMI^ 1 


111^ MMI^ 1 


111^ MMI^ r 




b ^ 












/ y 












/ / 

/ 


\ \ 
\ \ 


A"! 
// 1 








/ 

/ 


\ \ 
\ \ 


1 

1 








/ 


\ \ 
\ \ 


1 






I 




\ \ 












\ \ 
\ \ 


1 

ji 






Mid 1 


"J J II "J J 1 


\ Y 

III J J II II J J ii«" 


1 

J iiiiJ J 1 


iiJ J 1 





102 104 106 108 1010 1012 1014 1016 

e [eV] 

Fig. 7. a) Broadband electromagnetic emission produced at the source 
location via synchrotron and IC cooling by the electron distributions 
in Fig. [6J). The solid curve refers to the downstream component and 
the dashed curve to the upstream one. b) Expected flux at Earth for a 
source located at a distance D = 100 Mpc. The effect of photon-photon 
absorption during propagation has been taken into account. 



of 10^^ erg in the volume of the cluster The spectra are cal- 
culated under two assumptions: pure synchrotron losses (panel 
a), and synchrotron plus IC losses (panel b). The case in panel 
a is not realistic, but it allows us to understand the effects in- 
troduced by IC losses. Because of very hard injection spectrum 
of electrons below 100 TeV (see Fig.|4]l, the synchrotron losses 
lead to a standard E^^ type steady state spectrum of electrons 
with a high-energy cutoff which is determined by the cutoff in 
the injection spectrum. Although the production rates of pairs 
upstream and downstream are quite similar, more intense syn- 
chrotron losses downstream suppress the electron spectrum in 
that region. The effect of IC losses, which above 100 TeV takes 
place in the KN regime, leads to a more complex feature, in par- 
ticular to a hardening below the cut-off. Since IC losses affect 
the two sides of the shock in the same way, the slightly higher 
production rate of pairs downstream results in a higher flux of 
cooled electrons at energies up to 10'^ eV . At higher energies 
the synchrotron losses become dominant (see Fig. [Sj, and thus 
determine the spectral shape of electrons. 



4. Radiation Spectra 

We can now calculate the emitted radiation spectra based on the 
electron distributions in Fig|6]b. The broadband Spectral Energy 
Distribution (SED) of photons produced via synchrotron and IC 
emission is shown in Fig. |7^. Solid and dashed lines show the 
downstream and upstream contributions, respectively. Because 
of the enhanced magnetic field, the synchrotron emission is 
higher by a factor of * 10 downstream, as compared to up- 
stream. At the same time, the synchrotron peaks of the SED both 
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in downstream and upstream regions are located around 100 keV 
(slightly below and above 100 keV, respectively); the larger mag- 
netic field downstream is effectively compensated by the shift of 
the maximum of the electron distribution upstream to higher en- 
ergies (see Fig.|6][)). Due to the homogeneous distribution of the 
target photons for the IC scattering (2.7 K CMBR), the resulting 
IC gamma-ray spectra essentially mimic the electron distribu- 
tions. In particular, the Compton peak in the upstream region is 
noticeably shifted compared to the position of the Compton peak 
characterising the downstream region. 

Note that unlike previous calculations of proton spectra in 
galaxy clusters, we do not assume an a priori spectrum, e.g. of 
the "standard" form E^^exp[-E/E*], rather calculate its shape 
within the model of DSA, taking into account energy losses. 
Interestingly, our self-consistent calculation shows a less steep 
fall off in the final emitted radiation spectrum above the cut-off 
energy. The effect is the result of the combination of two factors: 
the accurate treatment of DSA correctly taking into account 
energy losses, //) the accurate calculation of the shape of the 
produced electron spectrum. In this regard we note that the com- 
monly used delta-function approximation significantly deviates 
from the results of accurate calculations jKelner & AharonianI 
I2008h . 

While some hints of non-thermal hard X-ray emission 
have been reported from Coma and Ophiuchus clusters (e.g. 
[pusco-Femiano et al 1999; Eckert at al. 2008), these objects 
have not yet been detected in gamma-rays. Our calculations 
show that if protons are effectively accelerated by accretion 
shocks, one should expect both hard X-rays and TeV gamma- 
rays from nearby and powerful clusters. The detectability of the 
X-ray and gamma-ray fluxes depend on the total energy of pro- 
tons accumulated in the cluster over the lifetime of the source 
(up to 10 Gyr) and the distance to the source. Also the angular 
size of the emission is likely to play a role, since nearby clusters 
of galaxies will appear as quite extended sources. Normalising 
the total energy of protons in the cluster to 10^^ erg, the lumi- 
nosities in both synchrotron and IC components of emission are 
expected at the level of L ~ lO'*^ erg/s. To estimate the expected 
flux, Fe = L/4-nd^, one has to fix the distance to the object. For 
example, in the case of the Coma cluster {d = 100 Mpc) the re- 
sults are shown in Fig. I?}), where the effect of gamma-ray photon 
absorption by the extragalactic background light (EBL) has been 
taken int o account. In the calcula tion we used the recent EBL 
flux from lFranceschini et al.ll2008l The energy flux for the syn- 
chrotron component appears to be of the order of 10 erg cm"^ 
s"', while the IC gamma-ray flux is approximately one order of 
magnitude less, given that gamma-rays with energy exceeding 
10 TeV are effectively absorbed in the EBL. The detection of 
these fluxes from extended regions such as clusters of galaxies 
is not an easy task, but possibly feasible with the next generation 
hard X-ray and TeV gamma-ray detectors. 
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Fig. 8. Proton spectra at the shock location for the case of a shock com- 
pression ratio R = 1. The solid curve refers to an acceleration time of 
10 Gyr, the dashed curve to an acceleration time of 5 Gyr The shock 
velocity is 2000 km/s, the magnetic field upstream is Bi = 0.3 yuG and 
downstream B2 = 1B[. 
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Fig. 9. Spectra integrated over the upstream (2 and 4) and downstream 
(1 and 3) regions for the proton distributions in Fig. [8] The solid lines 
correspond to an accelerator age of 10 Gyr and the dashed one to an age 
of 5 Gyr. 



5. Shock Modification 

We have shown in the previous sections that protons can be ef- 
ficiently accelerated at cluster accretion shocks. In this scenario, 
the reaction of the accelerated particles on the structure of the 
shock itself can be quite significant, leading to its modification. 
In turn, the modification of the shock reflects on the particle 
spectrum. This complex interplay leads to a variety of non-linear 
effects. A great deal of study has been dedicated to non-linear 
shock acceleration in recent years (we refer to Malkov & Drury 
I2OOII for a review). One of the consequences of shock mod- 
ification is the increase of the acceleration efficiency coupled 



with a significant hardening of the energy distribution of pro- 
tons. As a result, the available kinetic energy transferred to ac- 
celerated particles is accumulated at the highest energies close 
to the cut-off. This part of the spectrum provides the dominant 
contribution to the pair production, and ultimately to the broad- 
band electromagnetic radiation emitted by the secondary elec- 
trons. Therefore the features of the shock modification should 
be reflected in the spectrum of the non-thermal electromagnetic 
radiation. 

In order to explore this possibility, we take into account the 
modification of the shock in a simplified way, a full non-linear 
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Fig. 10. Pairs production spectra for an accelerator age of 5 Gyr. The 
solid line represents the downstream component, the dashed line refers 
to the upstream component. 
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Fig. 12. Electron spectra after cooling, downstream (solid line) and up- 
stream (dashed line), obtained by considering the energy losses in Fig. 
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Fig. 11. Energy loss rates for synchrotron plus IC losses. The magnetic 
field upstream is B] = 0.3 fiG (dashed line), in the downstream region 
= 7Bi (solid line). 



treatment being beyond the scope of the present study. We as- 
sume a compression ratio at the shock R - 1, which results 
from the linear theory of DSA when the relativistic component 
of the fluid (i.e. the accelerated particles) dominate the system's 
dynamics. Analogously, the magnetic compression factor is in- 
creased to ^ =7. In this case the spectral slope predicted by the 
theory in the absence of losses is a = 1.5. The calculated pro- 
ton spectrum at the shock surface, multiplied by E^'^ is shown 
in Fig. [8] for a source age of 5 and 10 Gyr. One can see in the 
spectrum, just before the cutoff, a bump which is induced by the 
energy losses due to interactions with the CMBR. 

In Fig.|9]we show the spatially integrated spectra of protons 
in the upstream and downstream regions. The curves are multi- 
plied by in order to emphasise the fact that the energy trans- 
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Fig. 13. a) Broadband radiation spectra produced at the source by the 
electron distributions in Fig. 1 12b. downstream (solid line) and upstream 
(dashed line), b) Energy flux at the observer location, after absorption 
in the EBL, for a source distance of 100 Mpc. 



ferred to protons is accumulated just before the cut-off energy 
region. 

The secondary pair production spectra are shown in Fig. [TO] 
while Figure [TT] shows the energy loss rates for electrons up- 
stream and downstream. Compared to Fig. |5] the synchrotron 
cooling in the downstream region is increased due to the higher 
magnetic field. This has an impact on the cooled spectra, as 
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shown in Fig. [12] The spectra are normalised to a total energy 
released in high energy protons of 10''^ erg. The spectral charac- 
teristics are similar to those shown in Fig.|6] but, due to the hard 
spectrum of protons, the amount of energy released in pairs is an 
order of magnitude larger than in the case of R - 4. 

This effect is directly reflected in the broadband radiation 
produced by secondary electrons (Fig.[T3b). The X-ray luminos- 
ity is increased by an order of magnitude compared to the case 
of R - 4. The same applies for the IC radiation emitted in the 
upstream region. 

Figure [T3b shows the expected photon fluxes at Earth, for 
a source located at a distance d - 100 Mpc. Because of the 
intergalactic absorption, the gamma-ray flux above 10 TeV is 
dramatically suppressed. This leads to the shift of the maximum 
of the arriving IC gamma-ray spectrum to 10 TeV with a flux at 
the level of 10"'^ erg/cm^s. 

6. Conclusions 

Proton acceleration in galaxy clusters was studied in the frame- 
work of DSA via a detailed time-dependent numerical calcula- 
tion that includes energy losses due to interactions of protons 
with photons of the CMBR. For realistic shock speeds of a few 
thousand km/s and a background magnetic field close to lyuG, 
the maximum energy achievable by protons is determined by the 
energy losses due to pair production and ranges from a few times 
10"* eV to a few times 10''^ eV. 

We performed the calculations assuming that acceleration 
takes place on time scales comparable to the age of the clus- 
ter Since steady state is never achieved in this scenario, a time- 
dependent treatment is required. Particle spectra, when calcu- 
lated including the effect of energy losses, exhibit interesting 
features. The decay of the spectrum above the cut-off" energy is 
not exponential. Its dependence on energy is shallower due to 
the flat profile of pair production timescales in the cut-off energy 
range. The time-dependent distributions of protons are used to 
calculate accurately the production rates of secondary electron- 
positron pairs. These electrons cool rapidly via synchrotron ra- 
diation and IC scattering which proceeds in the Klein-Nishina 
regime. The effect of the hardening induced by the KN cross- 
section is visible in the IC radiation spectra both in the upstream 
and downstream regions of the shock. For the fiducial Coma-like 
cluster used in this work, the synchrotron and IC peaks of the 
electron broadband SED are at comparable levels and the asso- 
ciated flux from a source at the distance of a 100 Mpc is expected 
at the level of 10 erg cm"^ s ' in the X-rays and an order of 
magnitude lower for TeV gamma-rays. Note that although the 
maximum of the gamma-ray emission is located above 100 TeV, 
it unfortunately cannot be observed due to severe intergalactic 
absorption. The expected gamma-ray flux from clusters of galax- 
ies is at the limit of the sensitivity of present generation instru- 
ments, however it may be detectable with the future generation 
of detectors. The optimum energy interval for gamma-ray detec- 
tion is between 1 and 10 TeV. 

The detectability of clusters in hard X-rays and gamma-rays 
associated with interactions of ultrahigh energy protons with the 
CMBR, depends on the value of the parameter A = W(,2/d^QQ, 
where W(,2 = W/10^^ erg is the total energy released in cosmic 
rays normalised to 10^^ erg, and t/ioo = d/lOO Mpc is the dis- 
tance normalised to 100 Mpc. Obviously, the best candidates for 
detection are nearby rich galaxy clusters like Coma and Perseus 
located at distances d ~ 100 Mpc. On the other hand, due to 
the large extension of the non-thermal emission (as large as sev- 
eral Mpc), and given that the angular size of the source cc l/d, 



the probability of detection reduces with the distance slower than 
l/d^. Nevertheless, as long as the total energy in accelerated pro- 
tons does not significantly exceed 10^^ erg, the visibility of clus- 
ters of galaxies in X-rays and gamma-rays is limited by objects 
located within a few 100 Mpc. 

The chances of detection of non-thermal emission of clus- 
ters related to ultrahigh energy protons, especially in the hard 
X-ray band, can be significantly higher if protons are acceler- 
ated by non-linear shocks modified by the pressure of relativis- 
tic particles. In this scenario, a large fraction of the energy of 
the shock is transferred to relativistic protons. Moreover, in this 
case a very hard spectrum of protons is formed, thus the main 
fraction of non-thermal energy is carried by the highest energy 
particles. These two factors can enhance the luminosity of X- 
and gamma-ray emission of secondary electrons by more than 
an order of magnitude, and thus increase the probability of de- 
tection of clusters located beyond 100 Mpc. 

In this paper we do not discuss the gamma-ray production re- 
lated to interactions of accelerated protons with the ambient gas 
which can compete with the inverse Compton radiation of pair 
produced electrons. The relative contributions of these two chan- 
nels depends on the density of the ambient gas and the spectral 
shape of accelerated protons. The flux of gamma-rays from pp 
interactions can be easily estimated based on the cooling time of 
protons, fpp X! 1.5 X 10'^/n_4 s, where «_4 = n/10"^ cm"-' is the 
density of the ambient hydrogen gas, normalised to lO"'* cm"^. 
Then, the energy flux of gamma-rays at 1 TeV is estimated as 
Fy(~ ITeV) « 6 X 10"'^A-W62n_4/dioo erg cm"^ s"', where k is 
the fraction of the total energy of accelerated protons in the en- 
ergy interval between 10 to 100 TeV (these protons are primarily 
responsible for production of gamma-rays of energy ~ 1 TeV). 
For a proton energy spectrum extending to 10'** eV, this fraction 
is of order of a- ~ 0.1. Thus for an average gas density in a clus- 
ter like Coma, n ~ 3 x 10""* cm"-', the gamma-ray flux at 1 TeV 
is expected at the level of 10"'^ erg cm"^ s"' which is compara- 
ble to the contribution of IC radiation of secondary electrons. In 
the case of harder spectra of protons accelerated by non-linear 
shocks, the contribution of gamma-rays from pp interactions is 
dramatically reduced and the contribution of secondary pairs to 
gamma-rays via IC scattering strongly dominates over gamma- 
rays from pp interactions. 
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